Electronic and optical properties of Cr and Cr-N doped anatase Ti02 
from screened Coulomb hybrid calculations 
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We studied the electronic and atomic structures of anatase TiCh codoped with Cr and N using 
hybrid density functional theory calculations. Nonlocal screened Hartree-Fock exchange energy is 
partially mixed with traditional semilocal exchange part. This not only heals the band gap underes- 
timation but also improves the description of anion/cation-driven impurity states and magnetization 
of the dopants. Cr and/or N doping modifies the valence and conduction band edges of Ti02 leading 
to significant band gap reduction. Hence, Cr, N and Cr-N doped Ti02 are promising for enhanced 
photoactivity. 

PACS numbers: 71.20.Nr, 71.55.-i, 61.72.Bb 
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I. INTRODUCTION 

For photocatalysis, titanium dioxide (Ti0 2 ) offers ex- 
cellent oxidation and charge transport properties to- 
gether with its wide availability, nontoxicity, and chemi- 
cal stability. For this reason, it finds many useful appli- 
cations such as photogeneration of hydrogen from water, 
degradation of pollutants under visible light irradiation 
and production of hydrocarbon fuels and dye sensitized 
solar cells (DSSC).[l|-[4[ Among the other polymorphs the 
anatase phase exhibits higher catalytic activity. [5| 

One important limitation is that anatase T1O2 has a 
wide band gap of ^3.2 cV [f| and can only absorb in 
the ultraviolet (UV) region (A<380 nm). This seriously 
reduces solar energy utilization to <~5%. Substitutional 
cation and/or anion modified titania has been proposed 
as an effective approach toget catalytic activity under 
visible light irradiation. |7Hl3l| Recent experiments have 
shown that Cr-N codoping in TiC-2 drastically enhances 
absorbance. dldH Moreover, Cr and Cr-N dopings were 
reported to increase the visible lig ht reactivity with in- 
creasing Cr incorporation, [lrl Il7| Ferromagnetism has 
also been observed in Cr doped anatase films. [Til. Il9| 
Theoretical studies predicted magnetization per Cr much 
larger than experimental findings. [13, Pure density 
functional theory (DFT) methods unreliably estimate 
half-metallic character for Cr doped Ti0 2 [l3, SI- For 
the optical properties, hybrid theoretical approaches has 
been recently used to get improved quan titative agree- 
ment with the experimental data. [l5l |23| Detailed elec- 
tronic structure investigations are still needed to get a 
proper description of band gap features for the doped 
systems. 

In this present work, we used screened Coulomb po- 
tential hybrid DFT calculations to investigate the modi- 
fications of the band gap properties of TiC>2 induced by 
Cr and Cr-N dopants and their effect on the correspond- 



ing absorption spectra. Formation energies have been 
calculated as a function of oxygen chemical potential to 
compare thermodynamical statbility of the doped struc- 
tures. We also obtained the magnetic moments and the 
charge states of Cr and N species. 



II. THEORETICAL METHOD 

For all systems, our calculations were performed based 
on the spin polarized hybrid density functional the- 
ory as implemented in the Vienna ab-initio simula- 
tion package (VASP).[25| Ionic cores and valence elec- 
trons were treated by projector-augmented waves (PAW) 
method. (2(| [27} The kinetic energy cutoff value was de- 
termined to be 400 eV. 

Pure DFT describes impurity and defect associated 
properties of transition metal oxides incorrectly. For ex- 
ample J Ti02 is predicted to be metallic in oxygen vacancy 
cases. [28( However, experimental studies show a semicon- 
ducting nature. This failure of DFT lies in the approx- 
imations to the exchange-correlation (XC) energy being 
usually local or semilocal. This flaw can be patched by 
using hybrid methods where a portion of the non-local 
exact exchange is admixed with the traditional semilo- 
cal exchange. We used the Heyd-Scuseria-Ernzerhof hy- 
brid XC functional (HSE06) [SHU] based on screened 
Coulomb potential reducing the self-interaction error. 
Technically, this results in a rapid spatial decay of HF ex- 
change which improves the convergence behavior of self- 
consistent iterations. The HSE exchange is derived from 
the PBEOjUjUl] exchange by range separation and then 
by cancellation of counteracting long range HF and PBE 
contributions, as, 
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where u is the range separation parameter for the screen- 
ing and a is the mixing coefficient. The parameter of 
w=0.2 A -1 is used for the exchange contributions as sug- 
gested for the HSE06 functional. |34| We determined the 
mixing ratio of the exact exchange to be 22% so that 
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FIG. 1: Relaxed geometries and total charge densities of substitutional N (a-c), Cr (d-f) and Cr-N (g-i) in anatase Ti02. 



calculated band gaps and lattice constants for both the 
anatase and the rutile polymorphs show good agreement 
with the experimental values. |13| 

In order to model the doped systems we constructed 
a 108-atom supercell by 3x3x1 replication of the con- 
ventional anatase unit cell. This structure is sufficiently 
large to accommodate spatial separation between the pe- 
riodic images of the defects. We traced all possible Cr 
and N doping configurations. In particular, For Cr/N 
codoping, the structure shown in Fig. [1] is found to be 
energetically the most favorable. For geometry optimiza- 
tion and density of states (DOS) calculations we used 8 
special fc-points to perform the Brillouin zone integra- 
tions. The calculated properties are converged such that 
they unnoticeably change when a higher fc-point sam- 
pling is used. We required a precision of 0.015 eV/A in 
the residual forces in every spatial component on all the 
atoms without fixing them to their bulk positions. 

The formation energies of the defects were calculated 



by using the formula, 

Ef = -Edopcd - -Epurc - n/iCr - «W + ™MTi + TO/Z , 

where Ed opc d (Epure) is the total energy of doped (pure) 
supercell. The chemical potentials of Cr, N, Ti and O 
are denoted by fici, /Un, MTi an d /io, respectively. De- 
pending on the presence of Cr (N) substitutional dopant, 
the number n (m) assumes values or 1. In thermody- 
namical equilibrium, the chemical potentials of Ti and O 
varies depending on the growth environment and must 
also satisfy the restriction /XTio 2 — MTi + 2^to ■ Under O- 
rich conditions, /io is taken from molecular oxygen, and 
then /iTi = /iTi0 2 — Eo 2 . Under Ti-rich conditions, ^Ti is 
derived from bulk (hep) Ti (/i T " lk ) and no is calculated 
from /xo = ^(A*Ti0 2 — A*Ti)- The remaining chemical po- 
tentials of Cr and N are taken from their natural phases 

(^Cr = ^{Ecr 2 3 ~ i^0 2 ) an d MN = |^N 2 )- 

Bader analysis quantitatively describes local charge de- 
pletion/accumulation. It involves integration of Bader 
volumes around atomic centers. These volumes are par- 
titions of the real space delimited by local zero-flux sur- 
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TABLE I: Average charge states (e) of dopants and their ad- 
jecent (nn) Ti and atoms from Bader analysis. 



Doping 


Cr 


N 


Ti(nn) 


O(nn) 


none 






+2.84 


-1.43 


N@0 




-1.38 


+2.31 




Cr@Ti 


+2.30 






-1.33 


Cr/N 


+2.29 


-1.13 


+2.81 


-1.19 



faces of charge density gradient vector field. We calcu- 
lated charge states of atomic species (in Table H]) using a 
grid based decomposition algorithm. [37| 



III. RESULTS & DISCUSSION 

The bulk properties of pure anatase were calculated 
using the 108-atom supercell. Our HSE-prcdictcd lattice 
parameters (a = 3.78 A, c = 9.45 A) and band gap 
value (3.23 eV in Fig. [2]) show remarkable agreement with 
the experimental data. Since nitrogen doping case has 
been discussed in detail elsewhere. [13f we will focus on 
Cr and Cr-N doped anatase to elucidate the role of Cr in 
photoelectrochemistry of Ti02. 

Cr- doped TiO-2 •' Experiments report incorporation 
of chromium at substitutional sites forming single- 
crystalline Cr-doped Ti0 2 .|||, H BS Bl They also con- 
clude that if chromium oxide forms it must be hi ghly 
dispersed and its size must be undetectably small. [1J- 
[Tsl HH In order to model the structure we substituted 
one Cr atom at a Ti site in the anatase supercell form- 
ing sixfold coordination with the nearest neighbor oxy- 
gens as shown in Fig. [T] All Cr-0 bonds are 1.90 A and 
slightly shortened relative to those of the undoped Ti02. 
The disturbance of the substitutional Cr on the lattice 
is very small agreeing with the X-ray diffraction (XRD) 
patterns. [l?], Hjl Under O-rich conditions, the formation 
energy of Cr@Ti doping has been calculated to be 1.06 eV 
(in Fig. [3]). These energetics are sensitive to the choice of 
the XC functional. For example, McDonnell et al. found 
a value of 0.36 eV with HSE06-PBEsol method. Un- 
der Ti-rich conditions the formation energy gets as large 
as 10.21 eV. Therefore, synthesis of substitutional Cr un- 
der O-rich conditions must be much easier as pointed out 
by previous calculations. [23f 

Substitutional Cr introduces 3d states which are dom- 
inant above the top of the valence band (VB) as pre- 
sented in the corresponding projected DOS (PDOS) of 
Fig. [2] For the minority spin component, fully occupied 
Cr 3d states appear isolated 1.07 eV above the valence 
band maximum (VBM). Our calculations are interest- 
ingly in agreement with the X-ray Photoemission Spectra 
(XPS) of Osterwalder et aZ. jT^] who reported the forma- 
tion of defect states at 1.0 eV above the VBM giving an 
intensity roughly proportional to the Cr concentration. 



HSE06 functional not only gives the relative positions 
of the defect states in the band gap but also better de- 
scribes the character of these states. For instance, un- 
like the full-potential linearized augmented plane wave 
(FLAPW) calculations of Ye et ai..(2fl| and LDA study of 
Peng et al., our HSE06 results predict Cr doped anatase 
to be semiconducting, not half-metallic. 

Chromium significantly modifies the conduction band 
minimum (CBM), too, by introducing empty Cr 3c? states 
at the bottom of the CB. Widely dispersing Cr- induced 
states effectively reduce the calculated band gap from 
3.23 eV to a value of 2.16(2.14) eV for the spin up(down) 
component. This band gap narrowing agrees well with 
the decrease of the binding energy of Cr 3d from 3.2 eV 
to 2.20 eV observed in the XPS spectra. 18 Depending 
of the choice of XC flavor, the band gap reduction was 
calculated by the previous theoretical studies to be 1.2 
eV,[H and 0.34 eVpU with FLAPW and HSE06-PBEsol 
methods, respectively. 

XPS core level and valence band data show that the 
majority of Cr is present in the 3+ charge state. [H, [TBI . 
[lH [Lsj In addition, the room temperature magnetiza- 
tion measurements reveal a ferromagnetic state for all Cr- 
doped anatase films, with a saturation magnetic moment 
of ~0.6 /xb per Cr atom.[l8[ Similarly, Zhang et al. re- 
ported Cr 3+ -associated ferromagnetism which increases 
up to ^0.42 ^e/Cr by lowering of oxygen pressures. fl9j 
Theoretical studies, on the other hand, predicted a charge 
state of 4+ with a magnetization of 2 • [III HU, [23| The 
discrepancy has been attributed to the possible presence 
of oxygen vacancies which leave excess charge in the sam- 
ples reducing the charge state from Cr 4+ to Cr 3+ . How- 
ever, without such a compensation, our HSE-calculated 
charge state of +2.3 and total magnetic moment of 1.21 
/l*b per Cr show a significant improvement over the ex- 
isting theoretical estimations. 

Cr/N-codoped TiO-2 •' Substitutional codoping of Cr 
and N in the anatase form has been confirmed by XRD 
measurements for various impurity concentrations. 
[I?], HH Initially, we considered all possible doping config- 
urations. After the relaxation, the Cr-N pairing shown in 
Fig-HIg-h) yields the lowest total energy. The lattice gets 
slightly disturbed because of the strong Cr-N interaction 
leading to a bond length of 1.62 A which is considerably 
shorter than the Ti-N bonding. The total charge density 
plots in Fig.[IJg) and Fig.[IJf) clearly shows the covalent 
character of Cr-N coupling while Ti-N bonding is more 
polarized. 

We calculated the formation energy for Cr/N codop- 
ing to be 4.43 eV under O-rich conditions. It rises up 
to 8.93 eV under O-poor growth environment as shown 
in Fig. [3] The experimental realization of Cr-N codop- 
ing is expected to be energetically viable under oxidizing 
atmosphere in agreement with previous studies. [i~5| 

The presence of substitutional N in the codoped 
anatase phase reduces the calculated magnetic moment 
to 1.00 fiB relative to the case of Cr monodoping. Be- 
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FIG. 2: Total (left) and projected (right) densities of states (DOS) of pure and doped (with Cr and/or N) anatase T1O2, 
calculated with HSE06 functional. Dashed (dotted) line indicates the Fermi energy (the VBM of pure anatase). 



sides, it insignificantly alters the charge state of Cr which 
is +2.29 as presented in Table. HI This indicates a strong 
derealization of N 2p states over TiC>2 bands. The UV- 
vis reflectance spectral measurements of Li et al. [HI and 
the X-ray absorption spectroscopy (XAS) of Chiodi et 
aZ.[l4| confirm the presence of Cr 3+ species. Our HSE- 
calculated charge state is reasonably closer than the pre- 
vious DFT+U predicted Mulliken charges on the substi- 
tutional Cr when N dopant is present. |23| 

Strong Cr-N pair interaction exhibits interesting fea- 
tures in the electronic properties in Fig. [5J First, the 
isolated impurity states associated with Cr 3d and N 2p 
in the (Cr@Ti) and (N@0) monodoping cases mutually 
passivate each other giving a clean band gap. This elim- 
inates the possiblility of photo-generated charge recom- 
bination due to well localized trap states. Secondly, the 
Cr and N concentration of ^0.93% in the computational 
cell causes a band gap reduction of 0.3 eV. McDonnell 
et al. calculated 2.66 eV band gap for Cr/N codoping 
by using HSE06-PBEsol XC functional. [3 Chiodi et al. 



reported the appearance of codoping-induced and Cr 3d 
dominant states at the top of the VB which translates 
to a band gap of 2.8 eV.[14| This observation nicely fits 
in our PDOS structure in Fig. [5] with an estimated band 
gap of 2.9 eV. For higher dopant concentrations, exper- 
iments report increased gap narrowing causing signifi- 
cantly red-shifted optical absorption. [14T - [l6j Thirdly, N 
2p states which appear in the monodoping case (N@0) 
in Fig. [2] significantly delocalize over the VB and the CB 
when Cr is present (Cr/N) due to strong Cr-N interac- 
tion. This PDOS feature shows remarkable agreement 
with the experimental observations, flil fl(| Therefore, N 
has a weaker effect on the optical response than Cr which 
dominantly modifies the band edges of Ti02- 

We calculated the absorption spectra of Cr- and Cr/N- 
doped cases in comparison with pure anatase (see Fig.[2|. 
Cr seems to be absorbing in the visible region signifi- 
cantly better than the Cr-N codoping. However, PDOS 
structure of Cr@Ti shows that the overall photocatalytic 
efficiency might be lower than the Cr/N case due to 
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FIG. 3: Calculated formation energies as a function of the 
oxygen chemical potential (relative to the value at its molec- 
ular gas phase) for Cr, N mono- and Cr/N co-doped Ti02 
structures.. 



the charge recombination in the presence of localized 
3d gap states. In fact, Cr-N codoped TiC>2 exhibits 
relatively higher catalytic reactivity under visible light 
irradiation. [16( 

IV. CONCLUSIONS 



The electronic, magnetic and optical properties of Cr- 
and Cr/N-doped anatase Ti02 have been investigated 
by means of screen Coulomb hybrid calculations. Hy- 
brid functionals incorporating exact exchange terms not 
only improve the description of defect states but also are 
promising to predict magnetization properties of d-band 
dopants. We have shown the dominant character of Cr in 
modifying the band edges of anatase. Substitutional Cr 
with or without N is very effective in band gap reduction 
causing significant red-shift of the absorption edge into 
the visible spectral region. Isolated N 2p states get de- 
localized over Ti02 bands due to mutual passivation by 
Cr 3d states giving a clean band gap for Cr/N-codoping. 
Thus, the enhancement of the optical absorption is emer- 
gent. Moreover, the calculations imply that the substitu- 
tional Cr and N pair reduces electron-hole recombination 
rate improving the overall photocatalytic performance. 
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FIG. 4: Calculated absorption spectra for pure, 
Cr/N-doped anatase Ti02 structures. 
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